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Abstract. We propose a simple model to study resonant tunneling through an 
organic molecule between to conducting leads, taking into account the vibrational 
modes of the molecule. We solve the model approximately analytically in the weak 
coupling limit and give explicit expressions for the thermopower and Seebeck 
coefficient. The behavior of these two quantities is studied as function of model 
parameters and temperature. For a certain regime of parameters a rather peculiar 
variation of the thermopower and Seebeck coefficient is observed. 

Although the model is very simple, we expect it to give some nontrivial insight 
into thermal transport properties through nan-devices. Furthermore, because we 
can provide an analytical solution, it may eventually serve as benchmark for more 
advanced analytical or computational methods. 
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1. Introduction 

In a recent experiment, Reddy et al. [I] studied transport through an organic molecule 
(benzenedithiol (BDT)) attached to a scanning-tunneling-microscop (STM) tip and 
a substrate, generated by a break junction using the STM tip. The authors were 
able to measure different transport properties, including the thermopower respectively 
Seebeck coefficient of this setup, which is is shown schematically in Fig. [TJ Quite 
generally, the idea to contact organic molecules to leads in a controlled fashion is rather 
fashionable and - if successful - can lead to a variety of interesting applications, ranging 
from nanoscale computing to quantum computing [2H1] . First progress of reproducibly 
and controlled contacting organic molecules has been made using carbon nanotubes 
[SHE], but similar attempts with smaller molecules usually suffer from several problems, 
among them the actual question of how the contact is established and how one can 
control its quality [SrtTT]. 

It is well known, that the physics of nanoobjects coupled to conducting leads is 
highly nontrivial (for an early review see for example [12j . more recent aspects are 
discussed in |13j). and already concentrating on electronic degrees of freedom only 
presents an extreme challenge to both experiment and especially theory. The latter is 
due to the fact that correlations on these nanostructures can usually not be ignored 
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Figure 1. Schematic experimental setup according to [Tj. The STM tip is 
connected to a reservior which maintains a constant tip temperature. The 
substrate, on the other hand, can be heated, thus creating a temperature gradient 
across the molecule. 



and hence one is plagued with complicated quantum impurity problems |14l I15j , which 
can be solved only numerically and only in equilibrium and for rather restricted model 
setups. The situation becomes even worse if one wants to take into account vibrational 
modes present in complex molecules [16 19 . Here, even advanced numerical tools 
rather quickly meet their limits [T5] , 

And last, but by no mean least, the conventional experimental probes used in 
nano-systems are transport properties, in particular by applying a bias voltage or, as 
in the above-mentioned experiment, a temperature difference across the nano-object. 
In this situation, the assumption of thermal equilibrium is quickly questionable, and 
reliable tools to study correlated systems out of thermal equilibrium - even only the 
stationary state - are not yet existing. 

Therefore, setting up simple models which preferably are solvable analytically and 
which allow to study the physical properties also off equilibrium in the stationary state, 
are still a highly relevant approach to obtain an at least qualitative feeling how various 
degrees of freedom influence the transport properties of nano-objects. Particularly 
thermal transport and here thermopower or Seebeck coefficients are an interesting 
area because possible applications in the area of energy conversion or storage have 
gained strong renewed interest over the past decade. However, in this area relatively 
little theoretical work on thermopower for mesoscopic systems exists to our knowledge 
so far HO] 

The situation in an experimental setup like the one shown in Fig. [T] has the 
advantage, that one can, for example by stretching the bond through removing the 
tip, control the coupling between the leads and the molecule very sensitively and 
actually shift the system into the weak coupling regime. Here, theoretical calculations 
become much easier. For our purpose, we will therefore assume that the system is 
in this regime, and suggest a very simple model which will allow us to see how the 
interplay between electronic and vibrational degrees of freedom influence transport 
properties. 

The model setup and general expressions for thermopower and possible situations 
will be presented in the next section. The central result, viz analytical expressions for 
the thermopower and Seebeck coefficient will be derivedin section [3j and a summary 
will conclude the paper. 
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2. Setup and model 



We study an organic molecule attached to two metallic leads. The leads are assumed to 
be in thermal equilibrium but can have different temperatures and chemical potentials. 
We will take the temperature and chemical potential of the right leads as reference 
and write Tr = T, Tl = T + AT respectively fJ,R = [J,, fiL = M + A/i. In the following, 
we use \i = as reference of energy. 

Using first-principle approaches, such a setup has been studied before in a quite 
general way |21j . concentrating however on the local properties of the molecule. 
Concentrating on a subset of degrees of freedom allows the use of more advanced 
analytical and numerical tools and has been applied to several model situations related 
to experiments mentioned in the introduction |18l I19| (see also further references in 
[19]). 

Here, we intend to investigate transport properties, and how local vibration modes 
influence them. Calculating transport for a general situation becomes rather quickly a 
cumbersome task. Therefore, in order to be able to gain some insight into the physics 
of the problem, while at the same time having as little parameters as possible, we 
concentrate her on a very simple model. Furthermore, restricting the study to a bare- 
bone model also allows for an explicit analytical solution. Such a solution is rather 
valuable in several respect, for example as possible benchmark for more advanced 
methods. 

To this end we assume that the HOMO of the molecule has a certain energy eo 
below the chemical potential of the leads, but the LUMO well above it, so that we 
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Figure 2. Schematic setup of the model. We consider only one relevant orbital 
of the BDT molecule, which we assume to be the HOMO. The coupling between 
the substrate respectively the STM tip are modeled as tunneling barriers. 

can neglect it. Also, all other electronic states of the molecule are assumed to be well 
separated, too. The model motivated above is shown schematically in Fig. [2j where 
we also explicitly included the various model parameters. 

Up to now the model is the standard setup to describe resonant transport through 
a mesoscopic system like molecules or quantum dots. Here, however, we also allow for 
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vibrational modes of the molecule, which we describe in a harmonic approximation, 
and which couple to the electronic degrees of freedom via a Holstein term. Although 
the calculation can be performed for several phonon modes, we concentrate here on 
the case of one single mode with frequency u>q. This is for example permissible when 
one studies low temperature, where only the lowest-lying modes play a role anyway. 
Finally, the molecule is coupled to the leads via tunneling barriers. 

Thus, the Hamiltonian for our model to describe the experimental setup from 
Fig. [T] reads 

H = Hl + Hm + #l-m (1) 

#L = J2 e «.k C a,k C a,k (2) 
k,a=L,R 

H M =e d}d + u> tfb + Aw d t d(6 t + b) (3) 



a=L,R 



(tJ a>B d + h.c) (4) 



where we use standard second quantized notation, A denotes the dimensionless 
coupling between the electronic and phononic modes, and t a the tunneling matrix 
elements to the leads. The operator c^ B denotes the projection of the band states to 
the barrier. 

Note that we did not include the spin degrees of freedom of the electrons in the 
model ([I]). Although this is at first glance a rather crude approximation, we want to 
employ it nevertheless for two reasons. First, including the spin makes already the 
molecular problem highly nontrivial and would prevent the semi-analytical solution 
we will discuss in the section about the results. Second, we are presently aiming at a 
more qualitative insight into the influence of vibrational modes on transport effects. 
To this end we believe that the additional complications due to additional electronic 
degrees of freedom will surely alter the results quantitatively, but possibly not strongly 
om a qualitative level. Nevertheless, calculations based on more realistic models are 
in progress. 

For our simple setup we can use the Meir-Wingreen formula [55] to calculate the 
electric current as 

J = Jo fdu [/ L H - f R (u>)] Pm (oj, A M , AT) (5) 



where pm(u, A/i, AT) denotes the electronic density of states on the molecule in the 
presence of a potential difference AV and a temperature difference AT between the 
leads, and 

/.(«) ' 
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the Fermi functions of the left and right lead, respectively. The quantity Jo finally 
collects the coupling parameters and the natural constants like charge and Planck's 
constant. 

Calculating the DOS in the presence of the leads for arbitrary coupling t a , A/i and 
AT is not possible. However, in the case t a — >• (weak coupling), one can approximate 
Pm(w, A/z, AT) by the result for the isolated molecule, which reads [2"3"] 

Pm (lj, A/i, AT) = -J W S ^ + A H) ( 6 ) 

I — — OO 
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and does not depend on A/i and AT explicitely. When we insert expression (|6]) into 
the Meir-Wingreen formula ([5]), the expression for the current becomes 

9 T 00 / 
6 i=-oo ' 

Since we are interested in the thermopower, the current has to be zero, i.e. for a given 
temperature difference AT between left and right lead we have to adjust A/x such that 

E W [/^o + a + AH)-Me + (Z + AH)]=0 (7) 
;=-oo ' 

In general, we are interested in low temperatures, i.e. typically we will have |eo|, 
ojq ^> ksT. In this case, only a certain range of values of I € [Zmin,'max], which have 
the property that eo + (/ + A)u; = 0(/cbT), will actually be important, because for 
the other terms the Fermi functions are either both zero or both one, i.e. cancel in the 
difference. This situation is schematically depicted in Fig. [3^. 

(a) (b) 
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Figure 3. Schematic representation of the different situations arising from 
Eq. (l7j. The left and right boxes denote the left and right leads, with the Fermi 
functions as full lines. In the left lead, we included some A/i and AT. In (a) 
there are several polaronic modes which lie in the Fermi window, while in (b) the 
temperature is much smaller than the level splitting, resulting in effectively only 
two modes within the Fermi window. 



Let us take this idea to the extreme and assume, that T is low enough such 
that only one vibrational mode with index Iq will contribute, see the schematic 
representation in Fig. [3]d. The necessary conditions for transport in this case are 

eo + (lo + A)w < 

< eo := e + (l + 1 + X)uj < uj 
from which can be determined. The equation ^ then simplifies to 

1 1 

-up — Afj. \ 



1 + exp 
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k B (T+AT) J 
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1 + exp 



\ k B T J 



1 



1 + exp 



( ep~A^ \ 
\k B (T+AT) ) 



1 + exp 



(JbJ) 



(8) 
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This equation can be easily solved for A/i for given model parameters and AT. 
Although this is strictly speaking not the case for the experimental setup in Fig. [T] we 
want to note that in nano devices eo can often be controlled through a gate voltage, 
€q = —eVa, and consequently the active polaronic mode respectively Iq is controllable 
exp er iment ally. 

3. Results 

Without vibrational modes, the thermoelectric response of the model has been 
studied before by Ermakov et al. 20J Within their model, the authors observed a linear 
relation between the temperature difference and thcrmovoltage. In the following we 
want to try to understand what modifications arise through the presence of polaronic 
modes in the molecule. 

In the present paper we want to concentrate on the situation shown in Fig. [3]3. 
Here we can use the fact that ujq 3> &bT and try to solve equation ^ approximately 
analytically and provide explicit expressions for A/i and the Seebeck coefficient as 
function of model parameters and temperature. 

The limiting cases eo — > and eo — > loq turn out to be uninteresting, because in 
the first case the combination/L(?o — w o) — /_r(?o — ^o) vanishes, while in the second 
/i(eo) — /.r(?o) —t 0. In both cases one thus recovers the solution by Ermakov et 
al.[2"0] in the limit when the level position is close to the chemical potential. 

The more interesting case therefore is eo finite and reasonably far away from 
the chemical potential. Then we can furthermore assume eo 3> k^T as well as 
| eo — w o| ^ ^bT, and replace the Fermi functions by f(x) f=s e~ l3x , if x > 0, respectively 
f{x) fa 1 — e l3x , if x < 0. Inserting these approximations into Eq. Q one obtains a 

Am 

quadratic equation for e'n (T+AT) . From the two solutions one must pick the one that 
has the correct limit 

AT 

Afi = (cj - e ) — 
as A —¥ 0. The solution fulfilling this requirement reads 

A ?£0^ 



A/i«(T + AT) In - 
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e k B (T+AT) 
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(9) 



Since this expression is rather hard to interpret, we inspect it behavior for small A. 
Taylor expanding it to first order results in 



. . AT A(T + AT) -n-2c r 
1 Iq + i 



gfc B T(T + AT) _ I 



The first term is identical to the result by Ermakov et al. [2D]. The correction is 
negative, i.e. the vibrations tend to initially reduce the thermopower. Note that the 
straight of this effect depends on the relation between eo and ujq: If ujq — 2?o < 0, it will 
be exponentially suppressed, while it is actually stronlgy enhanced for w — 2e > 0. 

Let us see in to what extent the approximate solution ^ is reliable. To this end 
we solve Eq. (18]) numerically for A/i. As generic parameter set we chose eo = wq/8 
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and X/(lo + 1) = 1/40. Due to the exponential structure a numerical solution tails for 
too small T. Using multi-precision arithmetic, one can extend the range of solvability 
slightly. We checked the stability of the root finding algorithm by comparing the 
results for A/x at A — > with the analytical result and found that stable solutions 
can be obtained for T > 5 • 10~ 2 wo using a precision of 100 digits. The resulting 
function A/z(AT) for T = 5- 10~ 2 u;o, 10 _1 o;o and wo/2 is shown in Fig. [i] as full lines, 
together with the approximate functions as dashed lines. For low T we observe a linear 
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Figure 4. A/Li vs. T _1 AT for various values of T. Other parameters see text. 
The dashed curves are the results from the approximate formula ([9l. 

dependence with negative coefficient, while at high T the sign of the thermopower 
reverses. For intermediate T, however, one can see a rather peculiar sign change 
of A/i as function of AT. Changing the model parameter does not influence this 
qualitative behavior as long as uiq — 2eo > 0. For example, decreasing ?o moves the 
zero in Afi towards smaller values of AT, but also decreases the magnitude of the 
effect. Increasing A produces a similar behavior. However, when uj — 2e < the 
features due to the vibrations vanish and one recovers more or less the behavior for 
A — > 0, i.e. linear with slope (eo — w )/T > 0. 

The curves resulting from the approximate formula give, as expected from the 
assumptions behind this analytical approximation, an accurate description for low 
temperatures, while for higher temperatures deviations become more pronounced. 
However, the qualitative behavior is well reproduced in all cases and we conclude that 
the analytical function represents a reasonable approximation for higher temperatures, 
too. 

Finally, one can look at the initial slope of A/i(AT — > 0), which, when divided by 
the electron charge, is the negative Seebeck coefficient S for the junction. The result 
as function of temperature for the same model parameters used in Fig. [4] is shown in 
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derived from the approximate formula ^ as dashed curve in Fig. [5] It apparently 
describes the overall behavior quite well, becoming asymptotically exact for low T. 
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Figure 5. Seebeck coefficient S = — QXtl as function of T for two 

characeristic values of en. Other parameters as in Fig. R] The dashed curves 
are the result from the approximate formula |10| . 

Interestingly it also provides the correct asymptotic for high T. Note that the unit of 
the thermopower is, as usual, determined by the universal ratio /ce/eo and is hence 
not model or parameter dependent. As already noted in the case of Afi, we here, too, 
obtain a sensible dependence of the Seebeck coefficient on en. When uj — 2e > 0, it 



shows strong variation with T due to the denominator in Eq. ( 10 ), while for u>o— 2en < 
the latter influence is exponentially suppressed and hence one obtains a temperature 
dependence which is very similar to the case without molecular vibrations. Since 
one can possibly control the value of eo and hence en experimentally, at least for 
electrically gated molecules, the described effects and their switching could in principle 
be observable. 



4. Summary 

A simple model for an organic molecule attached to two conducting leads via tunnel 
barriers was introduced. Main emphasis was laid on the inclusion of vibration modes 
for this molecule and to study their impact on transport properties, here especially 
the thermopower. We derived an analytical approximate formula for the thermopower 
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Afi(AT) and the Seebeck coefficient and showed that for a rather broad range of model 
parameters the influence of the molecular oscillations can strongly influence these two 
quantities. In particular, the thermopower can develop a pronounced extremum and 
sign changes as function of temperature gradient across the molecule, and the Seebeck 
coefficient showed pronounced nonmomotonic behavior in this regime. When one 
shifted the local electronic level of the molecule to larger values, one can observe a 
transition into a regime where the vibrational modes have negligible influence and 
the system behaves similar to the case without any elastic degrees of freedom. Since 
in some experimental setups the local parameters are accurately controllable via e.g. 
gate voltage or the like, we are certain that one can test these theoretical predictions 
to some extent. 

Of course, the simplifications introduced to allow for an analytical solution are 
rather severe: Electron spin as well as additional molecular levels were neglected, and 
the vibration modes were restricted to a single mode. Furthermore, for the calculation 
of the transport only two of the polaron levels were taken into account. The latter 
approximation can more easily be relaxed, even within the proposed model (see Fig. 
[3^,) . Similarly, including several vibration modes can be done straightforwardly (see for 
example 23 ), too, while taking into account more realistic electronic level structures 
will actually require more refined analytical approaces |16) or a numerical solution 
of the resulting model for the molecule [TU] ; and possibly also a reinspection of the 
.validity of the Meir-Wingreen formula for the current. 

In any case, adding these additional and physically important features to the 
model of the molecule makes the problem harder to treat. At least as long as we stay 
in the weak-coupling limit, we can expect that one can handle them reasonably within 
existing numerical approaches, or even, up to the solution of the transport equation, 
to a large extent analytically. The actual problematic part is hidden in the solution 
of the transport equation, even in weak coupling. Here, the exponential nature of the 
Fermi functions renders the numeric root searching extremely awkward and severely 
limits the accessible parameter and especially temperature regime (see e.g. discussion 
of Fig. Q. Finally, relaxing the condition of weak coupling to the leads leads to a 
full-scale many-body problem out of equilibrium, which at present cannot be really 
attacked with reliable analytical or numerical tools. 

Thus, to conclude, we believe that even such an extremely simple model will be 
of some importance to increase our understanding of properties of nano-systems, and 
can eventually also serve as a benchmark for more advanced methods to be developed 
in the future. 
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